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We discuss the role of neoclassical resistivity and local magnetic shear in the prediction of the 
Sawtooth period in tokamaks. When collisional detrapping of electrons is considered the value of the 
safety factor on axis, q(t, 0), evolves on a new time scale, r* = t v u* /(8y/e), where r„ = Ana 2 /[c 2 ?7(0)] 
is the resistive diffusion time, u f = u e /(e a ^ 2 ui t e) the electron collision frequency normalised to the 

O transit frequency and e = a/Ro the tokamak inverse aspect ratio. Such evolution is characterised 
■ 2 /3 

_ by the formation of a structure of size 5, ~ vj a around the magnetic axis, which can drive rapid 

evolution of the magnetic shear and decrease of q(t,0). In this paper we investigate two possible 
trigger mechanisms for a Sawtooth collapse corresponding to crossing of the linear threshold for 
the m = 1, n = 1 instability and non-linear triggering of this mode by a core resonant mode near 
the magnetic axis. The Sawtooth period in each case is determined by the time for the resistive 
evolution of the q profile to reach the relevant stability threshold. In both cases we focus on the 
s ! I scaling of the Sawtooth period as v* is varied. 
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I. INTRODUCTION 



When the safety factor, q, falls below unity, tokamaks experience a periodic oscillation in the plasma core in which 
core parameters exhibit a sawtooth-like waveform, with a relatively slow ramp-up of, for example, the electron tem- 
perature, followed by a rapid collapse. Understanding the characteristics of these Sawtooth oscillations is important 
for predicting the performance of ITER since they can degrade the core confinement, expel the alpha particles that 



c/2 , heat the burning plasma and couple to other instabilities that can severely limit operation. The length of the sawtooth 



period, which is terminated by the sawtooth collapse, plays a major part in determining the impact of the sawtooth on 



the tokamak performance. Following a sawtooth collapse, the radial profiles of the various plasma parameters evolve 



on transport timescales until a rapidly growing instability is triggered, producing magnetic reconnection on a fast 
timescale. It is this evolution, on the relevant transport time scale, that determines the sawtooth period. However 
the actual triggering event remains to be identified. In this paper we explore two potential avenues. 

The first, considered in Section 2, is in the spirit of the model proposed by Porcelli, Boucher and Rosenbluth [lj 
which, in particular, proposes a condition on the magnetic shear at the q = 1 surface for triggering the sawtooth 
collapse. This condition is loosely connected with diamagnetic stabilisation of the internal kink mode. Here, we 
establish more rigorously the stable window in an operating diagram for the m = 1, n = 1 drift-tearing mode and the 
. resistive internal kink mode, denned in terms of the plasma beta j3 = 8irp/B 2 (the ratio of kinetic and magnetic plasma 
\Q " pressure), and the instability drive, represented by a quantity A' which is inversely related to the potential energy for 
the internal kink mode SW. This diagram results from an earlier study Q of the stability of these two modes based on 
a plasma model with semi-collisional electrons and ions whose Larmor orbit exceeds the semi-collisional layer width 
where reconnection occurs. The evolution of the location of the core plasma state in this diagram, as the profiles 
7— I ■ of q and plasma pressure change during the sawtooth ramp, can be monitored and the triggering of the sawtooth 
crash identified as the point at which the linear stability threshold is crossed. This theory also predicts that the crash 
• rH ■ occurs when the magnetic shear, s = rq'/q, at q = 1 reaches a particular value, but is much more precisely defined 
than in the model of Porcelli et al. [lj. In particular, the electron temperature gradient, through the parameter 
r] e = dlog(T e ) / dlog(N e ) , is a significant parameter. Using a neoclassical conductivity the resistive evolution of the 
shear at q = 1 is monitored following a Kadomtsev reconnection of the q profile as a result of the sawtooth crash. 
The time to reach this critical shear then determines the sawtooth period; we deduce a scaling law for this. 

The second avenue, discussed in Section 3, builds on ideas of Park and Monticello [3] who realised that the resistive 
evolution of the q profile in neoclassical theory leads to a rapid, cusp-like drop of q on axis due to the effect of trapped 
electrons, with q rapidly falling to q ~ 0.8 in a typical MHD simulation. A more accurate treatment includes the 
collisional correction at small v* e = v e f e 3 ^ 2 u>te (here v e is the electron collision frequency, Ut e — vthe/Rq the electron 
transit frequency and e = a/R the inverse aspect ratio). This removes the trapped particle cusp behaviour very close 
to the axis so that dq/dr becomes zero there, although the value of q is still rapidly driven well below unity. In this 
case an alternative "trigger" for the Sawtooth collapse is taken to be q(0) < q cr it with q cr n ~ 0.75, or a similar value. 
The justification for such a simple criterion lies both in experimental data which tend to show, in many devices, such 
as TEXTOR, JET, and TFTR, for example, minimum values of qo ~ 0.7 just prior to a sawtooth collapse; and partly 
from the theoretical expectation that as the value of qo falls significantly below unity, new resonant MHD modes 



2 



become possible in this region of unfavourable average curvature: i.e. that the plasma core is increasingly at risk of 
undergoing a Taylor relaxation of the type suffered in Reversed Field Pinch devices [3]. We have investigated the 
consequences of this for the sawtooth period, based on the assumption that a core instability might occur when go 
falls to some critical value below 0.75, say, when the m — 3, , n — 4 ideal or tearing mode might be destabilised, for 
example. Such an unstable mode can couple toroidally to a mode resonant at q — 1, thus driving the m = 1, n = 1 
instability, which is invariably observed during the sawtooth collapse. Again the time for q to evolve to such an axial 
value provides a prescription for the sawtooth period for which we deduce again a scaling law. 

The scaling laws for the sawtooth periods resulting from these two models and their implications for ITER arc 
discussed in the Conclusions in Section 4 



II. LINEAR STABILITY AND THE SAWTOOTH TRIGGER 



In Ref. 0| we developed a unified theory of the drift-tearing mode and internal kink mode relevant to the m = 
1, n = 1 mode resonant at q = 1 in large hot tokamaks such as ITER. Specifically, we adopted a plasma model with 
semi-collisional electrons for which: 
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with = k y x/L s , k y the component of the perpendicular wavenumber lying within the magnetic surface, x the 
distance from the resonant surface and L s = Rq/s, the shear length. The resulting width of the electron current 
channel S is thus given by: 



6 = 



ky Vthe 



(2) 



We consider the ion Larmor orbit to be large, pi 3> S. The theory is characterised by two other key parameters 
: f3 = Q.h(3 e L 2 s / L 2 n and A' — s\/8W , where L^ 1 = —N~ 1 d r No e is the inverse of the equilibrium density gradient 
length, and A' is the instability drive @. For the m = 1, n = 1 modes, this is related to the potential energy of 
the internal kink mode, SW. The key results are that, at low values of /3, the drift-tearing mode, with frequency 
uj = oj* e (l + 0.73?7 e ), is stabilised by finite ion orbit and diamagnetic effects, provided that [see Eq. (42) of Q]: 



where 



with 



A' Pl < Ai, 



Ai = y/nfi log{A) - nfj(uj - 1)1, 



A = 



Cj = 1 + 0.73r? e , 



(3) 

(4) 

(5) 
(G) 



Co = uj/uj* e , u)* e = l/2k y VthePe/ L n , I is an integral defined in Q with approximate value ~ —r\ L J' 1 , and we have taken 
T e = Ti for simplicity. This result is valid at small A'pi ~ /3. At higher values of A'pi , i.e. as /3A'pi ~ 1 the unstable 
drift-tearing mode couples to a stable Kinetic Alfven Wave (KAW) until, at a critical value of the parameter f3A'pi, 
there is an exchange of stability, with the drift-tearing mode continuing at the same frequency but now stable, whereas 
the (previously stable) KAW becomes unstable. With continuing increase of $A' pi the KAW frequency drops towards 
Ld/uj*e — 1 and its growth rate also decreases, the mode eventually becoming stable when 



with 



A' Pi > A 2 



p l 2.7r] e l3 



(7) 



(8) 



where, Sq — y/oj* e i> e riL s /vthe, and r\ is the position at which q = 1. These results pertain to low /?, however, it has 
also been proved that the KAW is stable when (3 exceeds a critical value depending on r/ e , due to the effect of shielding 
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Figure 1: Stability boundary for the drift-tearing mode and kink mode as derived in Ref. Here So/pi = 0.01. 

of the resonant surface by plasma gradients. For the particular case, rj e = 2.53, this threshold is /3 = 0.34. Such a 
result is consistent with that of Drake et al.Q who found stability in the limit /3^(cj — 1) _1 ^>1. 

In Ref. [H, the stability limits for the dissipative internal kink mode were also determined. In particular, for 
SW < 0, one stability limit is given by, 



A' Pi < A 3 



where 



A 3 <5 = $ 
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0, hence 1\ = 1/2], while a general stability boundary for arbitrary 



[see Eq. (96) of Ref. 0, here we are using : 
ft, which ensured stability when f3 > /3 2 ~ \/5o/pi, was also determined. 

The combined effects of these stability boundaries is encapsulated in Fig. [1] In particular, when /3 -C 1, we observe 
the existence of two stable ranges for the stability index A'p^; namely 



A'pi < Ax, A!, < A' < A3. 



(11) 



Thus, as Fig. [T] clearly shows, we can identify a new window in f3 for the instability for the drift-tearing mode that 
was missed by previous calculations that exploited arbitrarily large j3. 

It is crucial to understand the relative magnitude of A' 2 and A 3 . In particular, for A' 2 /A 3 < 1, the low-/3 window 
of stability is accessible to the system. If we define/(?7 e ) = A' 2 /A 3 — 1, we can solve for f(f] e ) = 0, finding that the 
window of stability is present when f{rj e ) < 0. 

In Fig. [2] we show the plot of f(rj e ); we see that there is a critical electron temperature gradient for which 



A 2 /A 3 ^l, when Ve § V ^ . 

Some values of the new critical electron gradient r]^ are given in Table |T] 
From Eq. (96) of Ref. [2] , it is easy to show that 
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where s cr it is the critical shear at q = 1, and Q e the electron cyclotron frequency. This result can be compared with 
the result of Ref. [lj] , where a condition on the shear at the resonant surface is calculated from the inequality 
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Figure 2: Function / = f(rj e ) for ft 2 pi/5o = 10. 
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Table I: Critical electron temperature gradient ry e as a function of ft 2 pi/So, given by the solution of the equation 

f(Ve)=0. 



where is the linear growth rate of the mode and w« is the ion diamagnetic frequency. The main difference between 
Eq. (|13p and previous results is the explicit functional relation between critical shear and the energy 6W, and the 
precise electron temperature gradient dependence [20]. 

Once all the boundaries of marginal stability of these modes are determined, it remains to understand how they 
can be crossed. Now, noting the relationship Q 

A'ri = (15) 
6W v ' 



with SW = (rf / Rq)8Wt , where SWt is the energy calculated by Bussac et al. Q, §i is the shear at the q = 1 surface. 
The shear dependence of the key parameters ft and 5, reveals great sensitivity of the various thresholds to an evolving 
q[r) profile. In particular, assuming a parabolic pressure profile p — po(l — r 2 /a 2 ) and Tj = T e , (see Eq. (109) of Ref. 
[2|), for a given inverse aspect ratio e, we have 

h= . f° - 2 «0.5/K(n)] 2 , (16) 
e [aq'[ri)\ 

in JET or ITER. While it is desirable to run simulations of the Sawtooth cycle that couple stability criteria and 
transport evolution of all plasma prifiles, in this work we content ourselves with the analysis of the post-crash, 
neoclassical q evolution. 



III. EVOLUTION OF THE q(t, r) PROFILE DURING THE SAWTOOTH RAMP 

In this Section, we study the importance of the trapped particle correction to Spitzer resistivity in determining 
the duration of the Sawtooth ramp. During this period, which follows a collapse event, thermal equilibrium can be 
assumed to be rapidly re-established, but the current profile, and the q profile, evolve resistively towards a remote 
(and ideal MHD unstable) steady-state in which the toroidal current is = Eo/r)(r), with Eq constant, and rj(r) the 
resistivity. For the moment we ignore the effect of the Bootstrap current within Ohm's law. Neoclassical resistivity 
is given approximately by 0, [ltj| 

V(r) = Vs P (r)/(l - yfr/Rtf, (17) 

where -qsp is the Spitzer resistivity. With the electron temperature profile given by T e (r) = Tq(1 — r 2 /a 2 ) 4 / 3 , the 
Spitzer resistivity has the form 



( \ '/u 

Vs P (r) = — -J 

(1 — r l J a 1 ) 



(18) 
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We construct the relevant diffusion equation for the q profile in the cylindrical tokamak limit retaining one toroidal 
effect, namely the neoclassical correction to resistivity. Thus, 



dB e ,„ d , _ . d 

¥ = - C (VxE) 9=C -(,J z ) = - 



rjcf 



d 



Airr dr 



(rB e ) 



and using the definition of the safety factor, 



this becomes 
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dr \q 



ld_ 

r dr 



r) d 
r dr 



(19) 



(20) 



(21) 



where we have introduced the dimensionless variables defined by r = t/r r] , and r = r/a, with t v — Aira 2 /(rjc 2 ). The 
model for neoclassical resistivity is 



rj{r) 



(1 -r 2 )(l - v /ir 1 /2)2 ; 



(22) 



where e = a/Ro. Clearly, the fractional power in the trapped electron correction to Spitzer resistivity generates 
(unphysical) singular behaviour in Eq. pip, for r — > 0, i.e. in the vicinity of the magnetic axis. This is removed by 
including the transition from a neoclassical resistivity to Spitzer when 
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r 



3/2 



(23) 



(24) 



Incorporating this correction, the expression for the resistivity becomes 

V( r ) = 7=-a ' 

where = v e / '(e 3 / 2 u;t e ). In large tokamaks such as JET or ITER, the dimensionless parameter is extremely small, 
so that resistive evolution in the vicinity of the magnetic axis, though not singular there, is likely to be rapid: this will 
become evident from our numerical solution of Eq. (|21[) . Furthermore, although the scaling of the resistive diffusion 

time, t v (x a 2 Te^ 2 points to a much slower evolution of q(t, r) in ITER than in JET (possibly by a factor of ~ 60), the 
scaling of the small parameter, v* cx N e a/T 2 reduces this factor when considering core evolution times. For example, 
by expanding Eq. (f2"Tj) locally around x — 0, and employing Eq. (|2"4")l , one obtains the solution 



Qa(t) = g o (0)cxp (-t/T*) 



with 



R$N e 



'8V~e T} /2 



(25) 



(26) 



Hence, at early times, the safety factor undergoes an exponential decay on the timescale t*. Note that the presence 

2/3 

of the short timescale r* is the consequence of the formation in the q profile of a boundary layer of width 5* ~ vj a, 
which we assume to be destroyed by the crash itself, and thus not being present at t = 0, is developing only afterwards. 
We might suppose that a scenario with incomplete reconnection [ll|, cannot be affected by this time scale, lacking a 
mechanism to destroy the structure of width ~ v 2 J^ a. 



A. Evolution from a fully-reconnected Kadomtsev-like state 



As an initial q(r) state we first choose the fully reconnected Kadomtsev-like form. This is given by the analytical 
formula El 
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1 - tanh[(r - r mix )/<;} 1 + tanh [(r - r mix )/<;] 



<1k 



(27) 
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Figure 3: Kadomtsev-like pre-crash q in = qo/[l — r 2 /a 2 + (l/3)r 4 /a 4 ] and post-crash q pc profile as given by Eq. ([2"T)l . 



where qx — 1/(1 — 0.27r 2 ) is the Kadomtsev full-reconnected state which has been calculated numerically for 
initial qo = 0.75, qt n — Qo/[l — r 2 + (l/3)r 4 ] has been chosen for the the "pre-collapse" state, resulting in 
r mlx = y/9 - v / T44^o - 63/2 0.757, and <^ = 5 x 10~ 3 . 

The fact that we are using_a Kadomtsev-like post-crash state does not limit our conclusions. We could have chosen 
the Gimblett-Hastie state [12j, which differs from the one used in this work by having a Taylor-relaxed inner region, 
allowing for q < 1 around the axis. 

Since Eq. (|21[) is of the heat diffusion type, in order to solve it we need an initial value over the whole domain 
r 6 [0, 1], and two boundary conditions at r — and r = 1. The initial value is given by the function q pc defined 
in Eq. (|27|) . We then choose the Cauchy boundary condition q(t, 1) = qi n (l), i.e. the total plasma current is held 
constant, and the Neumann boundary condition <9 r (j -1 (t,0) = 0. The last one is chosen because we want our system 
to evolve, at very long times, towards an equilibrium magnetic field which is regular for r — > 0. To clarify this point, 
let us consider the case of constant resistivity. After setting dt = and integrating Eq. (|2Tj) twice, one obtains the 
equilibrium safety factor 

where C^a are two constants of integration. Then, if C2 ^ for r — > 0, we have q e q(r) — > C 2 _1 r 2 , which implies a 
divergent magnetic field, Bg/B = B ~ C^/r for r — > . Hence, we set C2 = 0. Then q eq (r) — > C^ 1 for r — > 0, which 
yields B/r ~ C\ for r — > 0. From this it follows that d r (B/r) — 0, or equivalently 9 r 5 _1 (i, 0) = 0. The result in Eq. 
(|28f can be easily generalised to the case of non-constant resistivity, giving 

edg \ l_ 



9e 9 (r) = 9pc (r = 1) ( / |^ ) TF ^. (29) 

'Ke) 

The integrals in Eq. ([29| can be performed analytically, then the result can be compared to the long time evolution 
of Eq. (|2"Tj) . Figure [7] shows this solution for i/r* = 120, 200, 300 and the analytical solution q eq (r) for 1/* = 10 -4 ; 
the analytical steady-state equilibrium is recovered. It is worth noticing that it is reached on a time which is much 
shorter than the resistive diffusion time, t v = 5 x 10 4 r*. Another important property of the steady-state solution 
of Eq. (|21| is that it must give a (radially) constant electric field [see Eq. (|T9|) ]. In Fig. [4b] we show the electric 
field E = fjr^ 1 d r (r 2 / q) calculated from the solutions in Fig. [3a] As expected, the final state is uniform through the 
domain of integration. Finally, in Fig. [5] we show an example of the evolution of the whole q profile, starting from 
the Kadomtsev post-crash state. The formation of the structure on the width 5* is evident. It is also clear that a 
very fast diffusion of the initial current sheet occurs at the mixing radius r m i X . In tokamaks, such a fast neoclassical 
evolution of the q profile seems to occur only when approaching stationary conditions Here we also note that 

the final equilibrium of Eq. ([2"9"]l is ideal MHD unstable to the m = 1, n = 2 mode, (since go < 1/2 @|) so, from an 
operational point of view, it should never be achieved! Incidentally, all the q profiles resulting from this analysis are 
very similar to those proposed within the incomplete reconnection model of Ref. [lj. 
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(a) The numerical solution of Eq. fl2 U for 
t/r» = 120, 200, 300. The analytical solution 
q eq (r) is approached from above as time 
increases. Here v t = 10~ 4 . 



(b) The electric field E = f)r~ 1 d r {r 2 /q) 
calculated from the numerical solution of Eq. 

(H) for t/r» = 120, 200, 300. The constant 
solution is approached as time increases. Here 
v, = 10" 4 . 



Figure 4: Long time evolution of Eq. (f2"Tj) for q profile and electric field. 



q(r/a,t/r, ) 
2 . 25 t 
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Figure 5: Solution of Eq. (J2TJ for time t/r, = 0, 2.5 , 37.5 , and 49.75 . Here u, = 0.001. 



B. Axial criterion 



In the previous section we drew attention to the fact that go undergoes a rapid downward evolution. One might 
consider what can limit such evolution of q on axis. It is known that, even before the ideal MHD m = 1, n = 2, 
instability becomes possible (when qo < 1/2), the tearing mode stability index A' m n for core resonant modes, such as 
m = 2, n — 3, or m — 3, n = 4, can become positive and potentially unstable, since diamagnetic stabilisation is likely 
to be extremely weak close to the magnetic axis, and average curvature is unfavourable [15]. Hence, it is tempting to 
look for a correlation between the onset of such modes, and the Sawtooth period. Since experimentally Sawteeth are 
observed to happen when qo ss 0.75, [16| we solve Eq. (|2ip for several values of ranging from = 10~ 4 to 2/* = .1 
and solve for the time at which q(0) = 0.75 . 
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0.01 
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0.09 


0.37 


0.005 


1.56 


0.075 


0.4 


0.0025 


2.43 


0.05 


0.49 


0.001 


4.75 


0.03 


0.62 


0.0005 


8.42 


0.025 


0.67 


0.00025 


15.6 


0.015 


0.09 


0.0001 


36.9 



Table II: Values plotted in Fig. H 




Figure [S] shows the results of Table [XT] graphically. A fit of the data gives 

tsaw fa = 7 x 10~ 4 + 0.034^ 2/3 + 0.004^ - 0.08^ 2 . (30) 
From this it is evident that when 0.001 < < 0.1, ts^^f/ 1 "* ~ > then 

tsaw ~ 7-,i/r 1/S oc Rl /3 N^ 3 T^ 6 sec, (31) 

where we do not distinguish between the length Rq and a, since JET and ITER share the same aspect ratio, Ro/a = 3. 
For smaller values of i/*, ts^vk scales as t*v~ \ thus the dependence in tsaw disappears. This is equivalent to 
saying that the fast evolution of q on-axis is a transient phenomenon regularising the q profile. For very small values 
of v* such that <5*/ri -C 1, the effect of v* on the global diffusion of the q profile is negligible. From a preliminary 
analysis, we also find that the presence of the Bootstrap current terms in Ohm's law reduces the strength of the 
electron trapping effect, i.e. the development of localised axial structures near axis. 

These results have to be compared to other similar ones. For example, Park and Monticello [3J, considering 

neoclassical diffusion within the q = 1 radius in nonlinear 3D MHD simulations, found tsaw oc T^ 2 Rq. Another 

interesting result is that of McGuire and Robinson pjj who found tsaw oc Ne^ 7 Te 9 ^ 14 without, or tsaw oc iVe^T 23 ^ 10 
with diamagnetic effects; these were obtained with an empirical fit to data. Notice that the weak temperature 
dependence we found in Eq. pip is mainly given by the different dependence of r\ and v* on temperature. In our 
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a = 1 m 


a = 3 m 


T e = 4keV 


Te = 25 keV 


t v ~ 400 sec 


Tr, ~ 24 x 10 3 sec 


r* ~ 0.86 sec 


r* ~ 3 sec 


v, ~ 0.01 


v,-6x 10" 4 


<5* ~ 4.6 cm 


<5* ~ 1.4 cm 



Table III: Resistive time t v , fast diffusive time t*, normalised electron collision frequency v t , and boundary layer <5* 

for both JET and ITER. 



case, when we consider the ratio t v /t*, we see that t^ et /t^ et ~ 10 2 , and t^ ter /t^ ter ~ 10 4 . If we take Ro/a = 3, 
and N e ~ 10 20 m -3 and compare the two machines, we obtain the results in Table [TTTT By using the results in Table 
Ull and Hill one obtains Tg E w 1.69 sec, and ~ 25 sec. A sawtooth period of 1.7 sec. is the longest observed in 

JET, while 40 sec. is the value empirically allowed to avoid triggering Neoclassical Tearing Modes |18j . 



C. Time evolution of shear at the q = 1 rational surface 



In Section 3, we stressed the fact that both f3 and A' , which are fundamental to the calculation of the boundary 
of marginal stability of the drift-tearing and kink modes, show a shear dependence. Hence, here we present the time 
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evolution of: n, the position of the q(r%,t) — 1 surface, the parameter (3 — 0.5/[aq' (ri)} 2 , the shear at the q = 1 
surface lj = riq'(ri), and finally, because of its relevance to core MHD stability, the value of q on axis. 
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Figure 7: The value of q on axis and the position of the 9 = 1 surface calculated from Eq. (|21[) for different v*. 
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Figure 8: The shear at the resonant surface s\ and the parameter j3 calculated from Eq. (p?Tj) for different 

From transport modelling simulations, we know that a typical value of critical shear at which the sawtooth is 
expected to be triggered is s± = 0.4 [19]. In our case, this is reached by the time the position of the rational surface 
r\ has approached the value at which g m = 1 in the initial Kadomstev pre-crash state. 



IV. DISCUSSION AND CONCLUSIONS 



In this work we have discussed the role of the dissipative, m = 1, n = 1, modes, in the prediction of the Sawtooth 
period in tokamaks, and explored the effect of neoclassical resistivity in determining the evolution of the plasma 
during the quiescent ramp phase of the Sawtooth. 
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In Ref. , we calculated the critical value of the stability index A[ l for crossing a linear stability threshold of the 
drift-tearing and dissipative kink modes, with gyrokinetic ions and semicollisional electrons. 

The stability thresholds derived in Ref. [H depend sensitively on the magnitude of certain plasma parameters at 
the (?=1 surface, such as the shear s\(t), and f3(r, t) = 0.5/3 e L^/L^. We have therefore explored the resistive evolution 
of these parameters. In addition, in order to address the possibility that the m = 1, n = 1 mode might actually be 
triggered by a core plasma instability near the magnetic axis, we have also monitored the evolution of qo(t). 

For the parameters defining the linear stability threshold of the m = 1, n = 1 mode, [i.e. ri(t), si(t), and (3(t),] 
we found negligible dependence on the collisionality parameter i>*, but faster evolution [in agreement with [3]] than 
would occur with Spitzer resistivity. On this basis, one would expect a scaling of the Sawtooth period from JET to 
ITER proportional to T^ 2 a 2 , assuming that energetic ion stabilising contributions are similar in both cases. 

However, if the rapid downward evolution of qo(t) were to be responsible for triggering a Sawtooth collapse, we find 
some sensitivity to the magnitude of z/* and a weaker scaling of the Sawtooth period with T e . These two different 
scalings may offer a means to distinguish the trigger mechanism in several machines. A suggested scaling from JET 
to ITER, in this scenario, is t$aw oc Rq^ 3 Ne Te ■ 
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